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ABSTRACT 

We discuss the potential of the eROSITA telescope on board the Spectrum-Roentgen- 
Gamma (SRG) observatory to detect stellar tidal disruption events (TDE) during its 4-year 
all-sky survey. We propose to distinguish such events from other X-ray transients using two 
criteria: i) large (more than a factor of 10) X-ray variation between two subsequent 6-month 
scans of the sky and ii) soft X-ray spectrum. The expected number of TDE candidates is 
~ 10 3 per scan (with most of the events being new discoveries in a given scan), so that a total 
of several thousand TDE candidates can be found during the 4-year survey. The actual number 
may significantly differ from this estimate, since it is based on just a few TDEs observed so 
far. The eROSITA all-sky survey is expected to be nearly equally sensitive to TDEs occurring 
near supermassive black holes (SMBH) of mass between ~ 10 6 and ~ 10 7 M o and will thus 
provide a unique census of quiescent SMBHs and associated nuclear stellar cusps in the local 
Universe (z < 0.15). Information on TDE candidates will be available within a day after their 
detection and localization by eROSITA, making possible follow-up observations that may 
reveal peculiar types of TDEs. 



1 INTRODUCTION 



It is widely accepted that supermassive black holes (SMBH), i.e. 
black holes with mass of hundreds of thousands to billions of solar 
masses, reside in the centers of most galaxies (see Ho 2008| for a 
review). However, apart from active galactic nuclei (AGN), where 
intense accretion of gas onto the SMBH takes place, and nuclei of 
nearby galaxies, where dynamics of stars and gas near the SMBH 
can be directly measured, such black holes remain hidden. A plau- 
sible scenario for revealing a compact object with M < 10 s Mjj 
in the center of an inactive galaxy is based on tidal disruption of 
stars passing sufficiently close to it (Hills 1975; Lidskii & Ozernoil 
1 19791) . According to analytical estimates, about half of the dis- 
rupted s tar's material should be cap t ured and th en accreted by the 
SMBH iGurzadian & Ozernoilll98ltlReeslll988l) . Numerical simu- 
lations confirm these suggestions predicting a flare of thermal radi- 
ation with t he peak at extreme ultraviolet (EUV)/soft X-ray w ave- 
lengths JAval, Livio, & Piranll200dlStrubbe & Ouataertll2009b and 
maximum luminosity of 10 43 - 10 45 erg s _1 . However, the expected 
rate of such tidal disruption events (TDE) is not high: depending on 
the stellar density in the nuclear cusp and the SMBH mass, it varies 
from 1(T 6 to 1(T 3 yr 1 per galaxy dWang & MenittI 12004). Con- 
sistent with these expectations, just a handful of TDE candidates 



1 SMBHs with M > 10 8 M o are expected to capture solar-type stars with- 
out disruption, since the corresponding tidal radius is smaller than the radius 
of the black hole's event horizon. Also, in the situation where the tidal and 
horizon radii are close to each other, general relativity effects cause the ca p- 
ture probability to strongly depend on the black hole spin (Kesden 2012). 



have been identified so far by X-ray (Komossa 2002; Donley et al. 
l2002l;lEsauej et al . 2008; Maksvm, Ulmer, & Eracleousll2010l) , UV 
JGezar i et al. 2009 and referen ces therein) and optical observa- 
tions ( Ivan Velzen et alj 1201 lbl and references therein). The ob- 
served spectral and temporal properties of these TDEs are consis- 
tent with expectations but the existing data are too sparse to provide 
stringent constraints on the theoretical models. 

Recently, it has also been recognized that powerful jets could 
emerge during accretion of stellar debris onto a SMBH, giv- 
ing r is e to radio and hard X-ray emission (Giannios & Metzge: 

|Le 



201 it Ivan Velzen. Kording. & Falckei r20l"H : iLei & Zhand |201 ll : 
Krolik & Piranll2012l) . The subsequent discovery by Swift of TDE 



candidates that fit well into this sc enario dLevan et al.l l201lt 
iBurrows et al.ll201ll ; ICenko et al.l ;2012) opened an extra dimension 
in TDE studies. It is currently unclear whether the observed rarity 
of such flares is owing to a high degree o f jet collimation o r to a 
low intrinsic probability of jet production dCenko et al.ll2012 ). Ob- 
viously, a larger sample of TDEs is needed to check the different 
possibilities and test the existing theoretical models. 

Deep wide-area surveys are especially powerful in discovering 
transient sources, usually by comparison with pre- or post- survey 
observations. As regards giant thermal flares from stellar TDEs, 
EUV/soft X-ray surveys are the most suitable. Up to several thou- 
sand TDEs had been ex pected to be detected d uring the ROSAT 
All-Sky Survey (RASS. ISembav & Westl fl993l) . but only a few 
were ac tually identified due to the lack of deep follow-up obser- 
vations JKomossa 2002; Do nley et alj2 002). Nonetheless, the con- 
tribution of RASS to the total number of TDE detections remains 
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Figure 1. Illustration of a TDE observation by eROSITA. The source falls 
into the telescope's field of view every 6 months, as indicated by the vertical 
lines. The average background level is shown by the blue dotted line. The 
green dash-dotted line marks the count rate that can be produced by Poisson 
fluctuations in the background with a probability P/, g =* 10~ 3 . The flux from 
a TDE candidate must exceed this level by a factor 5= 10 during at least one 
scan, which defines the TDE detection level shown by the red dashed line. 
See Appendix for further details. 



unrivaled. In the present paper, we investigate opportunities pro- 
vided by the upcoming Spectrum-X-Gamma (SRG) mission. This 
observatory is planned to be launched in 2014 to the L2 point of 
the Earth-Sun system and the first 4 years of the mission will be 
devoted to an all-sky survey in the 0.2-12 keV energ y band with 
the eR OSITA dMerloni et al]|2012l) and ART-XC JPavlinsky et"ZI 
[2012j) telescopes. The survey will consist of 8 repeated half-a-year 
scans of the entire sky, each with sensitivity several times better 
than RASS. Owing to this observational strategy, 'self-follow-up' 
observations of X-ray sources will be possible through scan-to-scan 
comparison. As one of the main TDE hallmarks is a giant ampli- 
tude of flux variations on the time-scale of years, searches for such 
events in the eROSITA All-Sky Survey (eRASS) data should be 
very efficient (see Fig.[]}. An additional identification tool will be 
provided by X-ray spectral analysis, taking advantage of the high 
sensitivity in the 0.2-2 keV energy range and the possibility of de- 
tection of harder photons with energies up to ~ 10 keV. Below, after 
a short overview of TDE X-ray emission properties, we specify the 
required identification criteria and estimate the expected TDE de- 
tection rate for eRASS. 



2 X-RAY SIGNATURES OF TIDAL DISRUPTION 
EVENTS 

The observational appearance of a given TDE may depend on the 
properties of the star (mass and type) and SMBH (mass and spin) 
as well as on the initial stellar orbit characteristics. To simplify the 
treatment, we shall use some 'average' TDE pattern that reflects the 
main observational properties of the currently available sample of 
X-ray selected TDEs. 

2.1 Temporal properties 

The flux decay of TDE flares is broadly consisten t with 
a power law with a slope of ~ -5/3 dKomossa 120021 ; 



IVaughan. Edelson. & Warwick||2004r) . This decay probably reflects 
the decreasing accretion rate, M, of the falling back stellar mate- 
rial by the SMBH, coupled with some constant radiative efficiency 
6 (the source luminosity L = eMc 1 ), since the viscous timescale 
is probably sho rter than the characteristic dynamical fallback time 
dPhinnevll989r) . However, this simple picture is probably valid only 
at a certain phase of the event, namely when M is below the critical 
value M Edd = £~ l L Edd lc 2 corresponding to the Eddington luminos- 
ity L Edd s= 1.3xl0 44 ( T $jjf-) ergs" 1 , but still above ~ 0MM Edd . At 

this stage, released gravitational energy is radi ated away by a geo- 
metri cally thin, optically thick accretion disk l lShakura & Sunvaevl 
1 19730 . For a solar type star with a peri-center distance R p = 3Rs 
(where Rs is the Schwarzschild radius), the peak mass rate occurs 
at time 

,5/2 
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M fl 



10 6 Mr. 



(1) 



after the instant of disruptio n tp and is highly supercri tical given 
6 = 0.1 and M BH < 10 7 M o dStrubbe & OuataertfcOOSl this delay 
is in fact highly sensitive to pericenter distance: t, oc R p ), There- 
fore, early on accretion takes place in radiatively inefficient regime, 
with an approximately constant luminosity L = L Edd emi tted by 
a thick accretion disk, which is likely subject to outflows dUlmen 
1999; Strubbe & Ouataert 2009). The boundary between the early 
(Eddington) and late (decay) phases lies at 
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yr 
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after disruption for R p = 3R S (r Edd oc R' . Istrubbe & OuataerJ 
2009). Alt hough the reality can be somewhat more complicated 
(see e.g. lLodato et al.[ 2009). we adopt the following approximate 
relation for the time profile of flares associated with TDEs: 
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L °(S) 
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for t < t 
for t Q < t < t { 

for t > t\ 



(3) 



where L qujes is the source luminosity in the quiescent state, Ln = 
V^Edd s> L quies (with t) being a geometrical dilution factor), and 
t\ = t + r Edd . Equation {3]l ignores the short period t, <k r Edd 

since Tjlr Edd =* 4x 10~ 4 (m s/7 /10 6 M o ) 21/1 ° (R p /3R s ) . However, 
as we further discuss in Section 1431 below, there is an interesting 
possibility to observe some TDEs during the rising phase of the 
flare, which lasts r,,, f < t, <k T Edd (under our assumption of short- 
ness of the viscous time compared to the fallback time) and can be 
significantly long for large M BH and/or R p . 

For SMBHs with M BH > fewx 10 7 Mq, only the closest stellar 
disruptions (small R p ) can lead to supercritical accretion, while the 
majority of TDE flares associated with such SMBHs are expected 
to be sub-E ddington (given a uniform distribution of pericenters, 
IUlmerMl999l) . Moreover, most of the emission will appear at UV 
rather than X-ray energies (see Section l2~2l below). We thus restrict 
our consideration to SMBHs with M BH < 10 7 Mq. 



2.2 Spectral properties 

An observer will detect only a fraction of the TDE bolometric flux, 
depending on the spectral energy distribution of TDE emission and 
the energy response of the detector used. Since this emission pre- 
sumably originates in an accretion disk around a SM BH, its spec- 
tral energy distribution can be approximated (IShakura & Sunvaevl 
1 19731) by a black body radiation spectrum with temperature 
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which corres ponds to the innermost regions of the disk (R ~ 5Rs , 
IUlmerlll9990 . Since this emission is fairly soft, the detection prob- 
ability must strongly depend on the absorption column density to 
the source, implying that X-ray selected samples might be biased 
against absorbed TDE events. No indication of significant intrin- 
sic absorption has been found from spectral analysis performed for 
known TDE candidates under the assumption of black body emis- 
sion. The estimated black body temperatures are KT\,b ^ 0.07 keV 
(the median value. iKomossal 12002 ; lEsquei et alj|2008i) . although 
the quality of the existing X-ray data hardly makes it possible 
to distingui sh a black body model from e.g. an absorbed power- 
law model JEsquei et al.ll2008r) . Power-law fits typically resulted in 
steep slopes (r > 3), i.e. the observed spectra are indeed soft. 

Disregarding the remaining uncertainty in the spectral shape 
of TDE emission, we use the (physically motivated) absorbed black 
body model in our calculations. This model has 3 parameters: ab- 
sorption column density N H , temperature KTbb and peak bolometric 
flux 



A 



4*4' 



(5) 



where d L is the luminosity distance to the source (see Section |4j. 
We freeze the first parameter at N H = 5 x 10 20 cirr 2 , which is close 
to the median value of Galactic absorp tion for random po sitions in 
the extragalactic sky as calculated bv lEsquei et al.l ( 120080 . In addi- 
tion to the Galactic absorption there might also be some intrinsic 
absorption due to edge-on orientation of the a ccretion disk and/or 
scree ning by unbound stellar debris material (IStrubbe & Ouataertl 
2009). However, in the former case the disk radiation might be di- 
rected preferentially away from the observer, and so the contribu- 
tion of such events to the net detection rate is likely to be small. In 
the latter case, the estimates of Strubbe & Ouataert (2009) indicate 
that unbound equatorial material builds up a 'vertical wall ' of de- 
bris with huge column density in the radial direction but subtending 
a small solid angle Afl ~ 0. 1 sr. Such screening will most likely re- 
sult in complete non-detection rather than detection of a TDE with 
a significant excess of absorption relative to the Galactic value. 

As follows from Eq. [4] the second parameter of the model 
obeys KT bb cc M B J . Unfortunately, there are only crude M BH esti- 
mates for known TDE candidates. The median observed temper- 
ature kTm — 0.07 keV suggests that typical black hole masses 
associated with TDEs are M BH « 10 6 M o . For comparison, for 
M BH = 5 x 10 6 and 10 7 M o one expects KT bb = 0.047 and 0.039 keV, 
respectively. We shall use these three fiducial values in our calcula- 
tions. 

As the accretion rate of falling back material decreases, the 
effective temperature of the radiation is also expected to decrease 
as K.T bb oc M' /4 , as long as the thin disk approximation remains 
valid ( IStrubbe & Ouataert 2 0091) . N onethele ss, late observat ions of 
five T DE candidates from RASS dVaughan. Ede lson. & Warwick 
2004) demonstrated no significant softening (and in fact indicated 
a marginal hardening for one of them) while the luminosities fell by 
more than two orders of magnitude. That may have resulted from 
transition to a different accretion mode accompanied by Comp- 
tonization of disk radiation in a hot corona and/or significant contri- 
bution from a relativistic jet. Taking into account these observations 
and because we expect eROSITA to detect TDEs mainly near their 
peak (see Section |4j, we disregard possible spectral evolution of 
TDE flares in our estimates below. 



It is also worth mentioning that somewhat higher temper- 
atures (i<Tt b ~ 0.12 keV) have been r eported for TDE can- 
didates detected in clusters of gal axies dCappelluti et al.l 120091 ; 
Maks vm. Ulmer. & Eracleousll2010f) . It is unclear whether this is 
owing to some selection bias, the specific physical environment 
and/or details of the spectral analysis. In any case, searches for 
TDE candidates with eROSITA in galaxy clusters will be chal- 
lenging due to the strong diffuse X-ray emission from the hot 
intracluster gas (which is much le ss of an issue for Chandra, 
iMaksvm. Ulmer. & Eracleousl 120101) . with a possible exception 
of brightest flares in nearby clusters (like the one reported by 
Ca ppelluti et alJl2009h . 



3 OBSERVATIONAL TECHNIQUE 

In formulating TDE identification criteria we should address the 
completeness of the resulting sample and reliability of the se lection 
proce dure. As usual in such cases (see e.g. Chapter 7 of IWall et al.l 
2003), if we wish to maximize the number of sources of the re- 
quired type, we should be ready to deal with significant contamina- 
tion of the sample by misidentified sources. In the case of eRASS, 
our primary goal is to obtain as large as possible a sample of TDE 
candidates that will make it possible to accurately measure the rate 
of such events in the local Universe. The associated systematic bias 
could then possibly be corrected by weighting the candidates ac- 
cording to their reliability. Hence, we choose rather loose identifi- 
cation criteria, based on key signatures of TDEs. 



3.1 eROSITA All-Sky Survey 

The SRG satellite will be rotating with a period of T = 4 hours 
around its axis pointed within a few degrees of the Sun, with 
the telescopes obse rving the sky at right angles to the axis 
I IMerloni et alj|2012f) . Following the orbit of the L2 point around 
the Sun the rotation axis will be moving at a speed of 1 deg per 
day. With the 1 deg-diameter field of view (FoV), eROSITA will 
complete an all-sky survey every 6 months and a total of 8 scans 
will be performed over the course of the mission. A typical posi- 
tion on the sky will receive ~ 240 seconds of exposure during one 
scan (and an order of magnitude more at the ecliptic poles). This 
exposure will be achieved by 6 consecutive passages of a source 
through the eROSITA FoV separated b y ~ 4-hour intervals (see 
iKhabibullin. Sazonov. & Sunvaevll2012l for further details). A typ- 
ical visit of a point source will thus last ~ 1 day, which is much 
shorter than the expected characteristic timescale of the flux decay 
during a TDE flare (few months). Hence, we shall regard such visits 
as single 240 s-long observations. 



3.2 Identification criteria 

For a TDE with a given spectral energy distribution, bolometric 
flux / (see Subsection l2.2t and redshift z, the count rate measured 
by eROSITA between 0.2 and 2 keV can be written as 

Ao.2-2 



Co.2-2(/,z) 



:/• 



(6) 



^0.2-2 (z) 

where Ao.2-2 (in units of counts cm 2 erg) is the flux-to-count conver- 
sion factor for a TDE at z — > and the correcting factor Xo.2-2 (z) 
accounts for the cosmological reddening of the spectrum. Both 
Ao.2-2 an d Xo.2-2 depend on the eROSITA response matrix and TDE 
spectrum. We have computed the redshift correction for different 
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Figure 2. Redshift spectral correction A"o.2-2 (z) for various spectral models. 



relevant spectral models using XSPEC dDorman & Ar naud 2001), 
and it proves to be significant (see Fig. [2}, especially for the model 
corresponding to M BH = 10 7 M Q . This has an obvious explanation: 
as the intrinsic TDE spectrum becomes softer with increasing black 
hole mass, a largely fraction of the photons red-shift out of the 
eROSITA energy band into softer bands (below 0.2 keV). 

Disregarding spectral evolution during the TDE, the count rate 
varies with time as Co2-i(f obsemr ) a L(t source (t obsmer )), where the 
relation t smlrcl ,{t obsen , er ) is governed by the source's redshift. 

Additional information on the nature of TDE candidates can 
be obtained by dividing the 0.2-2 keV energy band into several sub- 
bands, e.g. 0.2-0.4, 0.4—1 and 1-2 keV. The first one is sensitive to 
absorption along the line of sight, whereas the other two are well 
suited for estimating the effective spectral slope. We thus introduce 
a softness ratio SR = Cqa-i/Ci-2 to distinguish 'soft' and 'hard' 
sources. 

Taking into account the anticipated TDE properties (see Sec- 
tionO and eRASS characteristics (see Section [3TT1 . we propose to 
discriminate TDEs against other types of X-ray transients to be de- 
tected during the survey (such as AGN flares and GRB afterglows, 
see e.g. Komossa 2002) based on the following two signatures: 

• large (> 10) variation between the count rates (Co.2-2) mea- 
sured in two consequent half-a-year scans (see Fig.Q]; 

• soft spectrum (as indicated by the softness ratio SR), consis- 



tent with a weakly absorbed (Nr 
photon index r > 3. 



< 10 cm ) power law with a 



In addition, since we are dealing with soft X-ray flares, it is 
reasonable to limit TDE searches to Galactic latitudes \b\ > 30°, 
where interstellar absorption Nh < 10 21 cnT 2 . That will also min- 
imize contamination from cataclysmic variables an d flaring dM 
stars, whose appearance can resemble TDEs (see iDonlev et al.l 
120021) . 

Although the peak X-ray flux during a TDE can greatly ex- 
ceed the corresponding quiescent flux, i.e. the flux produced by the 
TDE's host galaxy, in reality the minimum flux will, in most cases, 
be dominated by the eROSITA background (see Appendix). There- 
fore, the first of our requirements formulated above implies that 
a candidate TDE must be at least 10 times as bright as the back- 



ground. Besides that, the count rate must be high enough for the 
softness ratio SR to indicate that T > 3 rather than T ~ 2 as typical 
of AGN. As demonstrated in Appendix, these conditions imply that 
at least Nim = 40 counts (i.e. 20 counts inside the half power diam- 
eter region) must be collected from a source between 0.2 and 2 keV 
during a 240 s exposure. The corresponding bolometric flux limits 
//,,„ for various spectral models are given in TableQ] Figure|3]shows 
a simulated eROSITA spectrum (for M BH = 10 s M G ) of a TDE with 
/ = fiim- Trial fits of this spectrum by various characteristic models 
demonstrate that it is feasible to reveal the spectral softness of a 
TDE using ~ 40 counts. 

We can finally formulate exact criteria that can be used for 
identification of TDEs during eRASS: 

1. Extragalactic (\b\ > 30°) location. 

2. At least 40 counts (or 20 counts inside the half power diame- 
ter region) detected from the source during the ~ 240 s exposure in 
a given scan of the sky. 

3. At least a factor of 10 higher count rate in comparison with 
the preceding or succeeding scan. 

4. Softness ratio SR > 2 in the bright state. 

The last criterion should not be used in application to relativis- 
tic TDEs, a subclass of TDEs discussed in Section|4~4l below. 



3.3 Detection in subsequent eRASS scans 

A TDE X-ray flare discovered in a given eRASS scan can remain 
detectable by eROSITA in its subsequent sky scans, i.e. 6 months, 
12 months etc. after the TDE trigger. This provides an interest- 
ing possibility of obtaining a crude (2 or more data points) long- 
term light curve of TDE emission. Specifically, to obtain a signif- 
icant (see Appendix for details) X-ray detection from a continu- 
ing TDE flare (whose position is known from its initial observa- 
tion by eROSITA), just ~ 4 counts during the ~ 240 s exposure 
are required. This implies, for instance, that a TDE flare will be 
detectable by eROSITA even if the X-ray flux has dropped by an 
order of magnitude during 6 months after the trigger. Such late de- 
tections, or even upper limits provided by non-detections, will be 
very useful for testing TDE models. 



4 PREDICTIONS 

The intrinsic rate, K, of tidal disruption events in the local Uni- 
verse, i.e. the number of TDEs per unit time and unit volume, 
should depend on a number of factors, such as the mass distribu- 
tion of SMBHs and the density profiles of nuclear stellar cusps, 
which are fairly uncertain at present. Hence, it is difficult to predict 
the TDE discovery rate for eRASS based on theoretical grounds. 
However, we can try to estimate this rate using the existing (poor) 
statistics of TDEs. Specifically, we can make use of an estimate 
based on the detection of two candidates i n the XM M-Newton 
Slew Survey and their non-detection in RASS dEsquei e t al. 2008): 
K ^ 5 x 10~ 6 yr~' Mpc~ 3 . This estimate is in line wi th other ob- 
servational results I IMaksvm. Ulm er, & Eracleous 2010) and with 
some theoretical expectations (Wang & Merritt 2004), although it 
may be affected by selection effects (see lGezari e t al. 2009). 
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Table 1. Summary of the results. The adopted TDE intrinsic rate ft is a crude estimate based on the currently available sample of X-ray selected TDEs (see 
text). It is assumed that all TDEs are associated with SMBHs of the same mass Mbh, for which three fiducial values are used. The detection limit is given in 
terms of the black body spectrum normalization / = Lo/(4ndi), where Lg is the bolometric luminosity and d^ is the luminosity distance to the source. The 
quantity N gives the total number of TDEs with the 0.2-2 keV flux above the appropriate detection limit. It is not corrected for those events that were already 
detected in the preceding scans (see the next column for the relative fraction of such events). The next two columns provide the relative fractions of events 
that are triggered in the supercrirical and rising phases of their light curves. The relative fraction of events for which eROSITA will provide at least 2(or 3) 
significant data points for constructing a long-term (> 6 months) X-ray light curve (see Section [331 is given in the last column. 
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Figure 3. Simulated eROSITA spectrum of a TDE at z - with flux near the 
eRASS detection threshold, i.e. providing 40 counts during 240 s exposure 
time. The adopted TDE spectral model corresponds to Mbh = 10 Af Q (see 
text) and the corresponding bolometric flux is 1.06 X 10 erg s cm- - . 
The spectrum is binned in energy so that the significance of data points is 
2cr. Also shown are best fits of the data by different models: solid - absorbed 
blackbody emission, dashed - absorbed power laws (Nh = 5 x 10 20 cm" 2 ) 
with the indicated slopes. An absorbed power law with T ^ 3 can be rejected 
with 99.5% significance based on^ 2 statistics. 



4.1 Detection rate 

We assume AC DM cosmolo gy with Q. M = 0.275, £2a = 0.725 and 
h = 0.7 ( K omatsu et alj|201 11 ). The luminosity distance as a func- 
tion of redshift is then given by 



d L (z) = a + z)d l 



Jo 



dz' 
E(z'Y 



(7) 



where d H = ■§- ~ 3.0/i ' Gpc 



4.3 Gpc is the Hubble distance, and 
H2 taogd [l99^. Defining J (z) = £ $., 



E(z)= V^,„(l+z) 3 +! 

equation (0 becomes d L (z) = d H (1 + z) J (z). 

Given the flux detection limit f Um corresponding to a given 
spectral model (Table [TJ, there exists the largest distance d L j im (or 
the corresponding redshift zu m ) from which a TDE flare at its peak 
could be identified with eROSITA: 



U> 



And L (zi im ) K(zum) 



: flu. 



(8) 



(1+Z«m) J(Zlim) K(ziim)- 



Aitd\f m 



(9) 



Solving this equation numerically with L = 1.3 x IO 44 erg s~' 
(the Eddington luminosity for a black hole with mass M B h = 
10 6 M o and geometrical dilution factor rj = 1), we obtain zn m = 
0.164. This is approximately two times further than the distance 
within which TDEs could be found with RASS (z ~ 0.09, 
iDonlev et all 120021 : lEsquej et al] l2008h . Similar calculations for 
M B h = 5x 10 6 and 10 7 M o , again assuming the Eddington luminos- 
ity at the peak phase, result in zn„, = 0.162 and 0.142, respectively. 
Remarkably, the limiting redshifts for the different SMBH masses 
prove to be similar. This happens because although the peak TDE 
bolometric luminosity increases with M B h, the spectrum becomes 
softer and most of the emission appears at frequencies below the 
eROSITA energy range. 

Consider now a small volume of the Universe at redshift z with 
radial width dz, covering an opening angle to on the observer's sky. 
The comoving volume of the slice is 

j if) 2 

dVAz) = codl-^-dz (10) 

E(z) 

(Hogg jl999l : see also Chapter 13 of |Peebleslll993h . Neglecting any 
evolutionary changes in the galaxy population between zu m and z = 
0, we can calculate the rate of TDEs in this volume as dN (z) = 
HdV c (z)dt. 

Taking into account both the maximum and decay phases of 
the light curve (eq. l[3}), a TDE flare remains detectable during a 
comoving time 

\3/5 



T lim (Z) - TEdd 



1 



\4nd 2 H f,, m (l+ z ) 2 J(z) 2 K(z)) 
Hence, at a given moment of observation there are 

dN(z) = RxT Um (z)xdV c (z) 



(11) 



(12) 



detectable TDEs. Substituting the corresponding expressions for 
dV c (z) and T/,„, (z) and integrating from z = to zn„,, we obtain 
the total number of detectable TDEs: 



N = a x I(zi im ), 
where 



3/5 



(13) 



(14) 



and 



I(Zl im ) 
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(zT 3l5 J(z) 4,s 



rutm K 

Jo Ti 



(l+z) 6/5 E(z) 



dz. 



(15) 



Finally, expressed in convenient units, the total number of de- 
tectable TDEs at a given moment is 

V-3/5 



"'^(ZtMstlw) 



Hi 



10- 



/(a,„).(16) 



4.2 Dependence on M BH and the phase of the TDE light curve 

Estimates of the S MBH mass fu nction 

JHopkins. Richards, & Hernquistl 120071 ; iGreene & Hoi |2007|) 
indicate that the spatial density of SMBHs is roughly con- 
stant in the mass range from ~ 10 6 to 10 7 M o but depends 
on redshift. The TDE rate may d emons trate some correlation 
(Brockamp, Baumgard t. & Kroupal 1201 lb or anti-correlation 
(Wang&Merritt 2004) with the central SMBH mass. There- 
fore, since the maximum TDE distance for eRASS is almost 
independent of SMBH mass, we may expect the TDEs found 
during eRASS to be distributed over the M BH range from ~ 10 6 to 
~ 10 7 M o . Since this distribution is currently unknown, we consider 
three different scenarios in which all TDEs detected during eRASS 
are characterized by a fixed M BH value: 10 6 , 5 x 10 6 or 10 7 M G . 
The resulting estimates are summarized in Table [TJ 

For M BH = 10 6 M G and adopting ft = 5 x 10~ 6 yr 1 Mpc~ 3 , 
u) = 2n (which corresponds to the \b\ > 30° sky) and R p = 3Rs , 
we find from equations ((2) and (Q~5} that tem = 0.1 yr and 



I(ZUm) 



0.0153, and finally from equation d!6t that 



N ~ 650 events can be discovered in the second eRASS scan 
through comparison with the first one. Figure |4]shows the distribu- 
tion of the TDE flares over redshift and as a function of time passed 
after stellar destruction. Approximately half of the TDE triggers 
are expected to take place during the supercritical accretion phase. 
For ~ 1/12 of the events, a second trigger (during the subsequent 
eRASS scan) will be possible, and in rare cases, three or more in- 
dependent f triggers will be possible. Thus, the majority of TDEs 
will be 'new' for each scan, i.e. not seen during the preceding one. 
For M BH = 5 x 10 6 and 10 7 M o , the fraction of 'old' TDEs, i.e. 
events observed in more than one eRASS scan is higher (Fig.[5]and 
Fig. |6j since the characteristic decay time positively correlates with 
the black hole mass (see Section^. 

Besides that, it turns out that the majority of TDEs discovered 
in a given eRASS scan will still be detectable (i.e. significantly ex- 
ceeding the background flux level) by eROSITA (see Section [3~3l l 
in the subsequent sky scan(s), i.e. 6 months, 12 months etc. after 
the initial trigger. Hence, eROSITA will typically provide 2-4 sig- 
nificant data points for constructing a long-term X-ray light curve 
of a TDE flare (see Tableland Figs.|4H6jl- As a result, eRASS will 
provide a unique sample of TDE light curves, which can be used to 
test TDE models. 

The above estimates were based on the assumption that the 
peak accretion rate exceeds the Eddington limit. In reality, only 
TDEs caused by sufficiently close stellar passages can provide a 
supercritical accretion rate, whereas more distant disruptions may 
be characterized by subcritical maximum accretion rates. Since the 
Eddington limit increases proportionally to the SMBH mass, the 
fraction of sub-Eddington TDEs may become significant for high 
M BH ; e.g. Ulmer (1999) estimate that the peak accretion rate ex- 
ceeds the Eddington limit in approximately 100% of events for 
M BH ~ 10 6 M o and in ~ 50% for M BH = 10 7 M o . In such a case, our 
estimates for the number of TDEs detected during eRASS should 
be reduced accordingly. 
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Figure 4. Distribution of TDE candidates over redshift for M B h = 10 6 M o . 
The total sample is shown by the solid black line. The events triggered as 
TDE candidates during the declining (post-supercritical) phase of the X-ray 
flare are shown by the two-dot-dashed black line. The events independently 
triggered as TDE candidates in more than one (two) eRASS scans are shown 
by the dot-dashed (dotted) black line. The blue, green and red lines denote 
events for which at least 2-, 3- and 4-point X-ray light curve can be ob- 
tained using data from consecutive eRASS scans. The thin vertical dashed 
line marks the redshift at which the numbers of events triggered during the 
supercritical and declining phases are equal. 



As was mentioned in Section [2Tl there is a non-negligible 
chance of detecting some TDEs during the rising phase of the X- 
ray flare, especially for large M BH and/or R p . We may place a rough 
upper limit on the number of such triggers: N rise IN E dd S T i7 T £rfrf> 
where Nem is the number of TDEs detected during the supercrit- 
ical phase. For M BH = 10 7 M o , up to ~ 20 rising TDE flares can 
appear in each eRASS scan (see Table [T). Since the rising phase 
is expected to last between a fraction of an hour and a few days, 
it might in fact be possible to see how the X-ray transient bright- 
ens up between its > 6 consecutive passages through the eROSITA 
FoV every ~ 4 hours. 



4.3 Revealing TDEs through eRASS-RASS comparison 

Apart from comparing the results of successive eRASS sky scans, 
TDEs can also be discovered by cross-correlating the data of 
the first eRASS scan with RASS dat^j. The characteristic un- 
absorbed 0.5-2 keV flux limit for RASS po int sources is 3 x 
10~ 13 erg s" 1 cnr 2 JBrandt & Hasingerl 12003) . Hence, to detect 
a factor of 10 increase in flux, a candidate TDE must reach a 
0.5-2 keV flux > 3 x 10~ 12 erg s" 1 cm" 2 during the first eRASS 
scan. This unabsorbed flux corresponds to a bolometric flux //,,„ = 
4.6 x 10~ 12 erg s~' ctrr 2 for the spectral model corresponding to 
M B h = 10 6 M o ). This limit is approximately 4 times higher than 
for TDE detection based on eRASS scan-to-scan comparison (see 



1 The sky cover age of XMM-Newton and Chandra surveys is too small for 
this purpose (see lBrandt & Ha singer 2005). 
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Figure 5. Same as Fig.g]but for M BH = 5 X 10 6 M o . 
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Figure 6. Same as Fig.|4]but for M B h = 10 7 M Q . 



Table [TJ. Solving equation (|9) for //,,„ results in zn,,, a 0.095. Us- 
ing equation l !16t . we find that cross-correlation of the first eRASS 
scan with RASS will provide an additional ~ 125 TDE candidates 
(assuming 7? = 5 x 10~ 6 yr 1 Mpc -3 , a> = In and M BH = 10 6 M o . 

In add ition, ~ 15 TDE candidates can be found using the 
method of iDonlev et al.l d2002f) . i.e. by looking for ROSAT sources 
with soft X-ray spectra that were at least > 10 times bri ghter during 
RASS than in subsequent observations. In their study dDonlev et al.l 
l2002h used ROSAT PSPC data, which cover only ~ 9% of the sky, 
for the follow-up, whereas with eRASS the analysis can be done 
for 50% of the sky (\b\ > 30°). Follow-up radio observations of 
such candidate s could be useful f or constraining models of TDE- 
associated jets JBower et al.l ( 120 12m see Subsection l5.2> . 



4.4 Tidal disruption events with jets 

Only two TDE fla res with signatures of a relativistic jet have been 
discovered so far dCenko et al.ll2012l) . Moreover, these detections 
were triggered in the hard X-ray (15-50 keV) band by the Burst 
Alert Telescope (BAT) on Swift. Therefore, it is difficult to predict 
how many such transients can be found during eRASS. 

A very rough estimate can be obtained by assuming that all 
TDEs with jets are similar to Sw J2058+05, the brightest Swift 
event of this type. Its estimated isotropic X-ray (0.3-10 keV) lu- 
minosity is L X j SO = 3 x 10 47 erg s -1 and the spectrum is con- 
sistent with a power law with r - 1.6 absorbed by a column 
density N H - 2.6 x 10 21 cm -2 , intrinsic to the source located at 
z = 1.185 l lCenko et al]|2012h . We simulated an eROSITA obser- 
vation of a source with such a spectrum at z = 1 and found that 
our count rate detection limit corresponds to an unabsorbed flux 
^0.3-10 = 2.3 x 10~ 12 erg s -1 cm 2 (we disregarded the K-correction 
since the spectral shape above 10 keV is unclear, whereas the spec- 
trum below 10 keV is variable and corresponds to K =s const for 
z i 1). Assuming a peak isotropic X-ray (0.3-10 keV) luminosity 
Lxjso ~ 5 x 10 47 erg s , the limiting redshift is zu m — 4.5. There- 
fore, TDEs with jets can potentially be found with eRASS even in 
the distant Universe. Clearly, the actual number of such events will 
strongly depend on the cosmic evolution of galaxies and their nu- 
clei. Ignoring the redshift dependence, the existing Swift statistics 
implies that roughly ~ 1 TDE with jets may be detected during 
each eRASS scan. 

An independent upper limit on the rate of TDEs with 
jets is provided by ROSAT data. The analysis of IDonlev et al.l 
(2002) revealed only one high-amplitude (> 20) flare with a 
hard spectrum (r < 2). This event is associated with the 
Seyfert 2 galaxy SBS 1620+545 dCarrasco et ail Il998l) at 7 = 
0.0516. On the WISE IR color-color diagram ( fwright et al.l 
|2010|) . the object falls on the boundary between Seyferts 
and spiral galaxies ([Wl-W2]=0.6, [W3-W2]=2.9, Data Tag: 
ADS/IRSA.Gator#2012/1017/03243 1.5650). The observed-frame 
absorption-corre cted 0.2-2.4 keV lu minosity is ~ 10 43 erg s _I in 
the bright state ( IDonlev et al.ll2002r) . which is much less than the 
observed luminosity of relativistic TDEs. Therefore, it seems more 
likely that this hard flare is an example of extreme variability in 
Seyfert 2 galaxies. Observation of SBS 1620+545 during eRASS 
will probably clarify the situation, since there have been no post- 
ROSAT X-ray obse rvations of the sour ce. Given that the same anal- 
ysis of RASS data ( IDonlev et alj2002l) revealed 5 soft flares - TDE 
candidates, we can tentatively conclude that the rate of relativistic 
TDEs is at least 5 times less than for 'ordinary' TDEs. Hence, per- 
haps not more than ~ 150 relativistic TDEs will be detected during 
each eRASS scan. 

In order to search for relativistic TDEs in eRASS data, one 
could use the same count rate criteria as for 'normal' TDEs (see 
Section l3"2l l but the spectral softness criterion should obviously be 
omitted or modified so that it selects hard rather than soft X-ray 
flares. 



1 This adopted value is somewhat higher than the luminosity measured for 
Sw J2058+05, but the corresponding SwiftfiQKT observations started some 
10 days after the beginning of the hard X-ray flare and so may have missed 
the brightest phase of the event. 
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5 DISCUSSION 

In order to reveal the true nature of the numerous TDE candidates 
to be found during eRASS, some additional identification/follow- 
up work should be done, as we discuss below. 



5.1 Cross-correlation with other surveys 

The localization accuracy for point sources detected during eRASS 
is determined by the eROSITA point spread function (PSF) av- 
eraged over its FoV (= 01°). For the faintest sources, it is ex- 
pected to be a ^ 12" (corresponding to a half-power diameter, 
HPD, of =a 29", lMerlonietal]|2012l) . whereas it will be much 
better for bright sources. At redshifts typical of our TDE sam- 
ple, z — 0.1, the corresponding linear size is thus < 24 kpc. 
Such accuracy should be suf ficient for searching for TDE host 
galaxies in optical (e.g. SDSS^, S PSS-Ill Collaboration et alj2012J. 
and Pan-ST ARRSl Kaiser et alj|201fjh and infrared (e.g. WISm, 
IWright et alj|2010 ) catalogs. However, association of a TDE can- 
didate with the nuclear region of a galaxy will be possible only for 
bright flares and/or nearby galaxies. 

Should the redshift of the likely host galaxy of a TDE can- 
didate be known, it may immediately indicate a very high lu- 
minosity of the transient and hence exclude its extranuclear ori- 
gin (i.e. a highly variable ultraluminous X-ray source). To ad- 
ditionally exclude s oft X-ray flares from active galactic nuclei, 
an optical spectrum dVeilleux & Osterbrocklll987l) and IR colors 
( IWright et al.ll2010l) of the counterpart galaxy could be used, if 
availably Therefore, cross-correlation efforts can greatly help in 
confirming/rejecting the TDE nature of candidates found during 
eRASS. 
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5.2 Follow-up observations 

The expected rate of 'TDE triggers' (see Section|4} during eRASS 
implies that there will be roughly two TDE alarms every day (data 
transfer from the spacecraft is planned to occur once a day). For 
some (perhaps a minority) of them, host galaxy identification will 
not be immediately clear. Optical follow-up observations for such 
events could clarify the situation and possibly also detect an op- 
tical signal associated with the TDE itself if carried out early 
enough after the TDE peak. Given the spectral energy distribution 
of TDE flares, it is more likely, however, to detect fading emis- 
sion in X-ray (e.g. with Chandra or XMM-Newton) and/or UV 
(GALEtf§ follow-up observations. Moreover, Chandra can provide 
excellent localization (< 1") and a high quality spectrum in the 
0.3- 7.0 keV energy band after a few kilo second observation (see, 
e.g.. lVaughan. Edelson. & Warwickll2004l) . An even higher signal- 
to-noise spectrum can be provided by XMM-Newton. A comple- 
mentary UV detection by GALEX would help constrai n TDE emis- 
sion m odels even better than with X-ray data alone I IGezari et al.l 
120091) . 

In addition, radio follow-up observations may be crucial for 
identification of relativistic TDEs, as has been the case for both 



1 http://www.sdss3.org/ 

2 http://pan-starrs.ifa.hawaii.edu/public/home.html 

3 Wide-field Infrared Survey Explorer, 
http://irsa.ipac.caltech.edu/Missions/wise.html 

4 Such flares, if found by eRASS, may be interesting themselves for study- 
ing AGN variability mechanisms. 
5 



On the contrary, searches for radio emission accompanyi ng 'or- 



dinary' TDEs have not yet given conclusive results (Bower et al.l 
l2012l : lvan Velzen et al.ll2012l) . 

We further note in this connection that blazar flares could in 
principle resemble relativistic TDEs. Indeed, blazar emission is 
also believed to originate in the vicinity of a SMBH from a rela- 
tivistic outflow directed at the observer. However, the two relativis- 
tic TDEs discovered by Swift were clearly different from known 
blazar flares and from any other known types of transients in that i) 
the X-ray luminosity was very high and ii) the broad-band spectral 
energy distributions (from radio to hard X-ray) were inconsistent 
with the so-called blazar sequence, in particular the optical to X- 
ray lu minosity ratio was very low (see Figs. 4—6 in lCenko et al.l 
120121) . Therefore, follow-up observations will likely be able to dis- 
criminate relativistic TDEs from blazar flares among hard X-ray 
flares detected by eROSITA. We also note that identification of 
new blazars in eRASS data is an important scientific task by itself, 
which can be addressed by different known methods , e.g. using in- 
frared colors provided by the WISE all-sky survey JMassaro et al.l 
l201ll) . 



Galaxy Evolution Explorer, http://www.galex.caltech.edu/ 



6 CONCLUSIONS 

We have demonstrated that eRASS will be very effective at finding 
stellar tidal disruption events. A uniquely large sample of several 
thousand TDE candidates can be obtained, making it possible to 
accurately measure the rate of stellar disruptions by SMBHs with 
masses from ~ 10 6 to ~ 10 7 M o and thus to obtain a census of 
dormant SMBHs and their associated nuclear stellar cusps in the 
local Universe (within z ~ 0. 15). Roughly half of the TDEs will be 
observed by eROSITA during the supercritical (Eddington) phase 
of accretion and the rest during the subsequent decay. In addition, 
there might be up to ~ 10 2 TDE flares caught during the rising 
phase, i.e. when the debris of the disrupted star have just reached 
the black hole. Furthermore, for the majority of the TDE candi- 
dates, eROSITA will still be detecting their fading X-ray emission 
in 1-3 subsequent sky scans, i.e. 6-18 months after the initial trig- 
ger. Hence, a unique sample of TDE light curves will be obtained, 
diminishing the need for dedicated X-ray follow-up efforts. There 
will thus be plenty of information on various stages of accretion 
of stellar material by SMBHs, setting stringent constraints on TDE 
theory. 

Importantly, TDE candidates can be recognized as such almost 
immediately (within a day) after their detection by eROSITA, so 
that prompt multiwavelength follow-up observations can be orga- 
nized. Such dedicated efforts may lead to the discovery of peculiar 
events among the many TDEs discovered by eRASS, e.g. associ- 
ated with a rapidly rotating SMBH or with the disruption of a rare 
type of star or a planet (such as the recently discovered hard X-ray 
transient IGR I12580+0134, lNikolaiuk & Waltej|2013l) . 

Of particular interest is the possibility of obtaining a sub- 
stantial sample of relativistic TDEs. Due to their extreme power, 
such transients can be detected by eRASS up to great distances 
(z ~ 4). The increased sample can shed light on the relation of 
jet-dominated TDE flares to 'normal' ones. Detection of relativis- 
tic TDEs at large redshifts will also provide a unique opportunity to 
study quiescent SMBHs (i.e. not in quasars) in the distant Universe. 
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APPENDIX 

Since we consider only inactive galaxies as possible hosts of TDEs, 
the eROSITA count rate associated with their quiescent (non- 
flaring) state will likely be below the background count rate, which 
is estimate d at 3.74 x IP" 3 c ts s" 1 arcmin" 2 in the 0.2-2 keV en- 
ergy range dMerloni et al. 2012JJ. This corresponds to 0.165 counts 
during a 240 s exposure inside the region of half-power diameter 
(HPD, 29") averaged over the FoV. Assuming Poisson distribution 
of the background counts, the probability for > 2 photons to be de- 
tected during this time inside HPD is P bg = 0.00066, i.e. < 10" 3 . 
Since we are looking for > 10 times amplitude flares above the 
background level, the detection limit for the source in the bright 
phase corresponds to 2 x 10 = 20 counts inside the HPD re- 
gion during 240 s. Therefore, the source count rate must exceed 
C/„„ = 0.167 cts s" 1 in the 0.2-2 keV energy band. The false rejec- 
tion probability of ~ 10" 3 is sufficient for our purposes, since the 
resulting sample of TDEs is not expected to exceed a few thousand 
events. 

Having 20 photons with energies from 0.2 to 2 keV inside 
the HPD region, it should also be possible to distinguish a soft X- 
ray spectrum as expected for a TDE flare from a harder spectrum 
typical of AGN. This distinction can be done in terms of a softness 
ratio,Sif = Co.4-i/Ci_2, the ratio of the count rates in theO.4-1 keV 
and 1-2 keV energy bands. This ratio depends on the absorption 
column density N H only weakly: for a power law spectrum with 
T = 2, SR(N H = 10 20 cm" 2 ) = 1.9 and SR(N H = 5 x 10 20 cnT 2 ) = 
1.6; for T = 3, SR(N H = 10 20 cm" 2 ) * 4.0 and SR(N H = 5 x 
10 20 cnr 2 ) as 3.3. Hence, having 20 counts within the HPD region 
and imposing the condition SR > 2, one can exclude at least ~ 50% 
of 'hard' (r = 2) sources but retain almost all 'soft' sources: for 
r = 3 and N H = 5 x 10 20 cnr 2 , SR> 2 with =* 80% probability. 
Since a typical TDE signal is believed to be even softer (r s 5 
in the 0.4-2 keV energy band), only a tiny fraction of true TDEs 
will be missed if the criterion SR > 2 is used. In fact, one could 
use an extraction region several times larger than the HPD region 
if the accurate shape of the PSF is known. Combining counts with 
appropriate weights could increase the signal-to-noise ratio by a 
factor of ~ V2 and also the efficiency of the softness ratio criterion 
(~ 90% of sources with r = 3 and N H = 5 x 10 20 cm" 2 will be 
retained). 

We also point out the possibility of exploiting the eROSITA 
sensitivity in the 2-10 keV energy range. For a source near the 
detection threshold, a 'soft' (r = 3 and N H = 5 x 10 20 cm" 2 ) 
spectrum will provide 0.30 cts in the HPD region in the 2-10 keV 
energy band for a 240 s exposure, whereas a 'hard' (r = 2 and 
Nh = 5 x 10 20 cm" 2 ) spectrum will provide 1.12 cts. The hard 
(2-10 keV) eROSITA background is estima ted at 0.04 cts in 240 s 
inside the HPD region iMerloni et alj|2012l) , although this number 
is rather uncertain due to the lack of data about the high energy par- 
ticle environment near the L2 point, which is believed to determine 
the eROSITA background in this energy range. Thus, a non-zero 
flux in the 2-10 keV range is expected for > 1/2 'hard' sources 
with a 0.2-2 keV count rate near the detection limit mentioned 



1 This estimate includes the contributions of the cosmic X-ray background 
(CXB) and particle background, with the former dominating below 2 keV. 
Below 0.5 keV, a major contributor to the CXB is diffuse background emis- 
sion from th e Milky Way who se intensity significantly (> 35%) varies 
over the sky (Lumb et al. 2002). This effect could be included in the real 
data analysis using t he ROSAT map of the Local Hot Bubble emission 
lSnowdenetalJl997l) . 



above. This additional information can be used to distinguish TDE 
flares from e.g. AGN variability. 

To summarize, we adopt a detection limit of 20 counts inside 
the HPD region during a 240 s exposure, i.e. C/,„, = 0.167 cts s" 1 
in the 0.2-2 keV energy band. We propose SR = C0.4-1/CV2 = 2 
as a boundary between 'soft' and 'hard' sources, which will allow 
one to retain as many as possible true candidates while rejecting 
a significant fraction of bogus ones. Of course, both the flux and 
softness ratio criteria can be adjusted if the background level or 
some other observational circumstances prove somewhat different 
from our assumptions above. 
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